Phonon line shape in disordered A^Cqo (A= K, Rb) 
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We present a calculation of the Hg phonon self-energy for a model of AaCso (A= K, Rb). The 
orientational disorder of the Ceo molecules is included, and the lowest order self-energy diagram 
is considered. The calculations illustrate that due to the violation of momentum conservation by 
the orientational disorder, Raman scattering can measure the decay of a phonon in an electron-hole 
pair, allowing the estimate of the electron-phonon coupling from such experiments. Comparison with 
experimental line shapes further provides support for a local correlation of the molecular orientations, 
with neighboring molecules preferentially having orientations differing by a 90° rotation. 



I. INTRODUCTION 



It is essential for the understanding of the alkali-doped 
fuUerenes AaCgo (A= K, Rb) to obtain reliable estimates 
of the strength A of the electron-phonon interaction. In 
metallic systems, such as AaCgo, a phonon can decay in 
an electron-hole pair. This decay leads to arLjextra width 
of the phonon. It was pointed out by AUenEJ that A can 
be estimated from the average of this extra width {"jph) 
over all q- vectors. This width can be measured in, e.g., 
neutron scattering.B It is, however, not easy to extract 
accurate values of A from such experiments for Cgo • Ra- 
man scattering provides an interesting alternative, since 
it gives a high resolution and since the Raman active 
modes {Ag and Hg) are the modes which couple to the 
electrons in the partly filled ti„ band. 

In the context of the fuUerenes, it was early pointed 
out that for the ordered system energy and momentum 
conservation forbids the decay of a q = intramolecu-. 
lar phonon in an intraband electron-hole pair excitation.0 
Due to the long wave length of the photons used in Ra- 
man scattering, this experiment would then not be ap- 
propriate for extracting A, since q « phonons are ex- 
cited. A3C60 has, however, strong orientational disorder, 
with each molecules taking-|essentially randomly one out 
of two preferred directions .□ Schliiter et al. therefore em- 
phasized that decay of a Hg phonon isjievertheless pos- 
sible, since q-conservation is violated,!, while the decay 
of a q = Ag phonon is not possiblefl It was later as- 
sumed that the violation of q-conservation is so efficient, 
that Raman scattering actually measures an average of 
Hg phonons at all q-vectors,Q as assumed in the Allen for- 
mula. Recently, however, the apnlicability of Allen's for- 
mula was anew put into questionB It was concluded that 
for a q— independent interaction, there is no broadening 
of the q = Hg mode, and that the broadening of this 
mode is entirely due to the q-dependence of the coupling. 
This led to a drastic reinterpretation of the coupling, con- 
stants derived from earlier Raman measurements.Ll Ex- 
plicit calculations have been performed for orientation- 



ally disordered systemsJjut the broadening of the q = 
modes was not studied.lij Therefore, although Allen's for- 
mula is the basis for deductions of A from Raman data, 
its validity for disordered A3C60 has to our knowledge 
never been explicitly tested. The purpose of this paper 
is to perform such a test. 

In Sec. II we present our model and the formalism, 
the results are presented and discussed in Sec. Ill and 
we give a summary in Sec. IV. 



II. MODEL AND FORMALISM 



We consider a model which includes the partly occu- 
pied tiu band. There are three iiu orbitals on each Ceo 
molecule i, which are connected by hopping matrix ele- 
ments t 
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The orientationaLdisordeiu has been built into the matrix 
elements iym,„/0 llj We want to describe the coupling 
to the intramolecular five-fold degenerate Hg Jahn- Teller 
modes. To describe the electron-phonon interaction, we 
use the Hamiltonian 
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where ujph is the a phonon frequency, bi^ annihilates a 
phonon with quantum number /i on site i, g is an overall 
coupling strength and V , are dimensionless coupling 



constantstjlla given by symmetry; 
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for pL — 1, .., 5 respectively, 
constant A is then given byl 
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electron-phonon coupling 
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where A^(0) is the density of states (DOS) per spin at the 
Fermi energy. We have assumed that the Hg modes are 
local Einstein modes with a local coupling to the elec- 
trons. Since these modes are intramolecular and since 
the interaction between the Cgo molecules is very weak, 
these assumptions should be very good. Actually, explicit 
calculations of the phonon dispersion curves for the Hg 
modes find that the dispersion is extremely small.E^I 

To describe the broadening of the phonon, we calcu- 
late the phonon self-energy due to the electron-phonon 
interaction. We consider the lowest order self-energy, ex- 
pressed in terms of the zeroth order electron and phonon. 
Green's functions. This is the formalism used by Allen.EJ 
In the present case, the phonon energy is not much 
smaller than the bandwidth and the validity of Migdal's 
theorem is therefore questionable. Contributions beyond 
the lowest order self-energy diagram may therefore be 
important. The on-site Coulomb interaction is large in 
these systems, which are believed to be close to a Mott 
transition.!!!! For simplicity, we have here nevertheless ne- 
glected corrections to Migdal's theorem and effects of the 
Coulomb interaction. 

We diagonalize the electronic part H'^^ of the Hamil- 
tonian. The eigenstates are labelled by n, the eigenval- 
ues are e„ and the corresponding zeroth order electron 
Green's function is G%{lo). The phonon self-energy is 
then 

n.,o'(^) ^~^Y.Y.l '^9Zra9iaG'i{u^')Gl{uj+0j'), 



(5) 



where a = (i/i) labels a combined site and phonon degen- 
eracy index and g"„j is the matrix element of the coupling 
to the phonon a between the one-particle states n and 
m. The frequency integral can be performed analytically, 
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where /„ = fi^n) is the Fermi function, i] an infinitesi- 
mal positive number. In the present calculation, we used 



77 — 0.0004 eV, comparable to experimental resolutions 
and much smaller than Uph and jph- 

The phonon Green's function D is given by 



D-^ = [z?"]-i-n. 



(7) 



where a matrix notation is understood and the noninter- 
acting phonon Green's function is given by 
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We introduce the reduced Green's function in momentum 
space 



where N is the number of molecules in the cluster, c!f j, 
the scattering coefficients with proper symmetry factors 
corresponding to the specific experiment, as discussed in 
the Appendix. 

We compute the coefficients c!f j, for Raman scattering 
experiments. We consider electronic transitions within 
the tiu and tig bands, assuming that the photon energy 
is larger than the ti^-tig splitting. To obtain a simple re- 
sult, we make the approximation that the hopping matrix 
elements are small. Calculations using realistic hoppings 
give qualitatively similar results. Details of the derivation 
are given in the Appendix. Then, we have an expression 
for cjf^ which only depends on fi, v and molecular orien- 
tations oil) at sites i, j, (o(i) = ±1). 

Between molecules with the same orientation, i.e., 
o(*) = o(j). 



4^^ = 1 (all m), 
For o(i) = l,o(j) = -1, 
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and for o(i) — —1, 0(7) — 1, 
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Finally, 



c!f„ = 0, otherwise. 
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In the following we are interested in D^,^(q = 
0,0;), since in Raman scattering the coupling is to 

Eu -0/^,^(^ = 0'^)' and in 



i?'-H-) = ^EE^M,.(q,-)- 
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In the q = phonon Green's function momentum con- 
servation is considered, to the extent that it survives the 



orientational disorder. On the other hand, the usage 
of ij'o'^^' instead impUes an assumption that the orien- 
tational disorder has completely destroyed momentum 
conservation and that we therefore can average over all 
momenta as in neutron scattering. The spectral func- 
tions are defined as the imaginary part of the Green's 
functions, i.e, p^(q, lj) = -(l/7r)ImL)^,^(q, lj). 

Allen's formula for the phonon broadening "fph (full 
width at half maximum (FWHM)), can be derived from 
Eq. (^ in the limit of Wph <C W and by ignoring the 
momentum conservation, and it is given by 
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The phonon level shift estimated from the real part of 
the self-energy can be expressed as 
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for a flat DOS. The second term in the bracket is the 
correction of order {LUph/W)"^. 
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FIG. 1. q=0 and local phonon spectral functions at differ- 
ent rotational energies. The local spectral functions (dashed 
lines) show structureless profile with the weak dependence on 
the rotational energy. The g = spectral functions signifi- 
cantly depend on the disorder structure of the fee Ceo lattice. 



III. RESULTS 

We have performed calculations for clusters with 256 
atoms on an fee lattice and with periodic boundary con- 
ditions. We take into account the fact that neighbor- 
ing molecules tend to anti-alignjlj due to the -Lonping 
integral being stronger for such configurations O'tj We 
control the randomness of orientation by introducing 
an Ising-type nearest neighbor antiferromagnetic inter- 
action for the two preferred orientations, i.e., Erot = 
J^~^^^i jyo{i)o{j) with Af the number of nearest 
neighbor pairs. First we start from a truly random ori- 
entation and then anneal the system to any desired ran- 
domness by using a ficticious temperature. Since the fee 
lattice is frustrated, there are many "antiferromagnetic" 
(AFM) states which minimize Erot at Erot.mm = —1/3. 
Therefore even a perfectly annealed system retains a cer- 
tain disorder, as can be seen in the electron density of 
states. Then the phonon spectral functions are averaged 
over many sets of orientational configurations. 

The phonon spectral functions for the scattering vector 
q=0 and the local (summed over all q) spectral function 
are plotted in Fig. H for different orientational disorders. 
The phonon line widths for q=0 (solid lines) and local 
(dashed lines) spectral functions are about the same over 
a large range of disorder (different values of Erot)- This 
shows that the orientational disorder effectively breaks 
momentum conservation, and it supports the assump- 
tion that the Raman scattering experiment can give an 
estimate of the electron-phonon coupling. On the other 
hand, a calculation for a periodic "ferro-magnetically" 
ordered system gave two unbroadened peaks due to en- 
ergy and momentum conservation, as will be discussed 
later. 



The validity of Allen's formula, Eq. (|15D, is explic- 
itly tested by changing the phonon frequency Wph- Since 
the self-energy is trivially proportional to g^, we concen- 
trated on the linear dependence of "fph on Wph for a fixed 
g. The numerical estimates for the broadening ^ph and 
the level shift Aoj are computed by fitting the q=0 spec- 
tral functions to a single-pole phonon spectral function 
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where A is the total spectral weight Jla Although the nu- 
merical spectral functions are far from being single-pole 
structured, the peak position and the FWHM, jph, can 
be understood in an average sense. 

The widths and shifts obtained from these fits are com- 
pared with the results of Allen's formula in Table 1. In 
Allen's formula we have used iV(0) = 6.2eV~"'^, which 
is the appropriate value for our model of A^Ceo. The 
coupling constant g was set to 0.03 eV. We find that 
Allen's formula is reasonably well reproduced from the 
q=0 spectral functions. The discrepancies between the 
numerical and Allen's estimated values range within 20- 
40 %, which we attribute mainly to finite size effects. The 
phonon broadening decreased as the number of sites was 
increased. Also note that jph shows a poorer agreement 
for small values of tUph, where finite size effects become 
more pronounced due to the small number of particle- 
hole pairs with excitation energy less than ujph. Due to 
computational limitations, we have considered a maxi- 
mum number of 256 sites although the convergence with 
respect to the system size is not good enough for a sys- 
tematic finite size scaling analysis. 



We have calculated the Raman scattering spec- 
tra for systems close to the local minima of anti- 
ferromagnetically ordered orientations by using realistic 
parameters for the coupling constaat g obtained from 
the photoemission spectra analysia22l. We then discuss 
the relevanoa of our calculation to Raman scattering 
experiments^ We have performed the calculations for the 
first two Hg-modea at frequencies of 0.0336 and 0.0542 
eV (270, 432 cm~^, respectively) and have taken aver- 
age over 10 sets of rotational configurations for Erot — 
-0.323, close to the AFM minima at Erot = -0.333. 

For the first two Hg modes the fitted 7pft's are 13, 47 
cm"-'^ for u}ph=270, 432 cm~^, respectively. These val- 
ues differ by almost a factor of 2 from the experimental 
estimates at 20 and 21 cm~^, respectivelyJj This illus- 
trates the fact that photoemission and Raman scattering 
give rather different estimates of the distribution of the 
coupling strength for the low-lying phonons.llZI The calcu- 
lated Raman spectrum for Hg{2) mode is shown in Fig. ^. 
The total spectrum can be resolved in two major peaks 
from the first two (// = 1,2 at A) and the last threecHl 
(/i ~ 3, 4, 5 at B) modes of the phonon. This leads to 
a spectrum which agrees reasonably well with the low- 
energy part of the spectrum in Ref. O, in the sense that 
the spectrum is strongly skewed towards higher energies 
and has adjiitiopal structures. This provides support for 
the beliefO'llSlla that there is a local 'antiferromagnetic' 
correlation, since the agreement with experiment is bet- 
ter than for the random orientation. The two narrow 
high frequency peaks in the experimental spectrum seem, 
however, to be beyond our treatment. 

To discuss the phonon level splitting into a doublet 
(/i = 1,2) and a triplet (^ = 3,4,5), we define the re- 
duced phonon self-energy at q=0 in analogy to the re- 
duced phonon Green's function i?^,^/(q = 0,a;) as 
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which is the corresponding self-energy for D^^i (q — 
0,0;). (See Appendix.) The real part of 11^,^(0, Wp/i) 
agrees very well with the phonon level shift previously ob- 
tained by fitting the spectral function for £'^,^(0, uj). The 
level shifts are better understood by investigating the 
imaginary part of n^_^(0,w), which, related to the real 
part by the Kramers-Kronig relation, gives the excitation 
spectrum responsible for the perturbational level shift. 
Spectra of fmn^_^(0, uj) extend from to the bandwidth 
and show different shapes for the doublet and triplet. 
For orientational configurations close to the AFM min- 
ima, those for the doublet are skewed to higher energy 
giving the smaller contribution to the level shift while 
those for triplets are skewed to lower energy leading to 
the larger level shift. For random orientations, the two 
spectra show essentially the same shape as the system 
size is increased. 

The imaginary part of 11^^^(0,0;) represents the spec- 
tral distribution of coupling of /i-th phonon to electron- 



hole pair excitations. More specifically, 11^^^ (0,0;) can 
be written as the energy spectrum of the state for the 
uniformly distorted field, ji/)^^''), due to /i-th Jahn- Teller 
phonon mode; 



lmn^,^(0,c.) = ^^ (a,/3|0(^)) 5{u-eo.+ep), 
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where the summation is over the electron-hole pair states 
\a, (3) {Sa > ep, £f3 < ^F with the Fermi energy ep) and 



i</''^^>= E '^L.kII'^^vIo), 



(20) 
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with the non- interacting Fermi sea |0). ipima is the anni- 
hilation operator for the orbital pointing along the m-t]\ 
principal axis, not along the molecular orientation, i.e., 

^iraa = V'jmcr for o{i) = 1 and f/ii^cr = Ipiycr^i^iya = -Ipixcr 

and Tpiza- = ^iza for o{i) == — 1. The energy expectation 
value of the distorted field can be related to the first mo- 
ment of lmn^.^(0, oj) as 



(0(^)|0(A')) 



= Eo 



/q°° a;lmn^,^(0, uj)duj 
/Q°°lmn^,p(0,o;)dw 



(21) 



where Eq is the ground state energy of non-interacting 
electrons. As expected near the AFM local minima, 
(0('^)|7J"='|0(^))/(0(^)|0('^)) has smaller value for the dou- 
blet and larger value for the triplet. 

The Jahn- Teller phonon couples to the distortion of the 
electronic density. For the doublet modes, the phonons 
couple to the density distortion along the principal axes, 
e.g., for /i = 1, {h\^ + bn) couples to (n„ + riiy - 2niz), 
the distortion elongated or squeezed along z-axis. Sim- 
ilarly, the triplet phonons couple to distortions along 
diagonal directions in xy,yz, zx-planes e.g, for /i = 3, 



(&I3 + bis) couples to {uis: 



where x,y axes are ro- 



tated from x,y-axes by 45° along z-axis. Therefore, the 
distortion for /x = 3, 4, 5 tend to distribute more elec- 
trons along the in-plane diagonal directions, which are 
the strongest bonding directions for ii„-orbitals in the 
fee latticeJl3 Furthermore, upon the inter-site coherence 
close to AFM minima, the distortions from neighbors add 
up constructively along the strongest bonding direction 
to give the energy gain for the triplet phonons, leading 
to the skewedness toward the low energy. 

It has been suggestedH that the locally broken cubic 
symmetry, due to the orientational disorder, could lead 
to a splitting of all five levels in the Hg modes. However 
within our model, it seems unlikely that the local energy 
splittings survive after averaging over all sites. Rather, 
the calculation indicates that the local energy splittings 
contribute to the broadening of levels, not to collective 
q=0 modes which survive disorder. We observe that the 
first two peaks in the experiment of Winter et al. could 



reflect states in isolated Cgo, since their observed energy 
shifts from the undoped Ceo values are too small com- 
pared to estimates of the phonon self-energy from the 
lowest order diagram and sometimes even have different 
signs. Considering the narrowness of these peaks, they 
could result from a many-electron molecular states, due 
to effects neglected an this study, and with a very weak 
coupling to the surrounding. 
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FIG. 2. Phonon spectral functions for realistic parameter 
Wph = 0.0542 eV and g = 0.0359 eV for Hg{2) mode. The 
thick solid line shows the q=0 phonon spectral function as can 
be observed in Raman scattering. The spectrum has underly- 
ing structures which can be resolved to two peaks A and B by 
plotting projected spectral functions for each phonon modes 
(thin lines). The local phonon spectral function (dashed line) 
shows no such underlying structure. For comparison, the peak 
positions for a completely ordered periodic lattice are marked 
by arrows. 

Our results are in strong contrast to the results of 
Aksenov and Kabanov.El They studied the problem in a 
model where the orientational disorder was described by 
a relaxation rate, leading to an imaginary self-energy in 
the electron Green's function. However, the breaking of 
the translational symmetry was not explicitly included, 
and momentum conservation was assumed in evaluat- 
ing various diagrams. The orbital degeneracy was also 
neglected. It was concluded that for a q-independent 
electron-phonon coupling there is no broadening of a Hg 
phonon and that the broadening is entirely due to the 
dependence of the coupling on the angle of q. This is in 
strong contrast to the results above, where we find that 
a q-independent coupling gives a substantial broadening, 
provided that the violation of the translational symme- 
try is included explicitly. As was discussed in Sec. II, the 
q-dependence of the coupling should be very weak, and 
the effect discussed in Ref. Ig not very important. 

Aksenov and KabanovS were not able to obtain the 
angular dependence of the coupling constants and there- 
fore they could not deduce the absolute magnitude of 



the coupling constants from Raman experiments, but in- 
stead estimated the relative coupling constants from Ra- 
man experiments. To obtain absolute estimates, they, 
used coupling constants deducedcJ from photoemissioncll 
for a free Cgo molecule and arbitrarily assumed that the 
coupling constant to the lowest mode is correct. Had 
they instead, equally arbitrarily, assumed that the cou- 
pling to the third mode is correct, all values of A^ would 
have been about a factor jUf two larger and if they had 
assumed that the deducedE2l coupling to the eighth mode 
is correct, all couplings would have been zero. 



IV. SUMMARY 

We have demonstrated that because the orientational 
disorder breaks the momentum conservation, Raman 
measurements of the Hg modes in the AsCgo should see a 
broadening due to the decay in electron-hole pairs. This 
enables Raman scattering experiments to derive electron- 
phonon coupling strengths. Furthermore, the explicit 
calculation shows that the width of this q=0 mode is not 
very different from the q-averaged width, where the latter 
is described by Allen's formula. We further find that the 
local " anti- ferromagnetic" correlation of the molecular 
orientation is important for obtaining the non-symmetric 
line-shape seen experimentally. This provides support for 
such a local correlation actually taking place in the real 
system. This leads to a large splitting between a two-fold 
and a three- fold degenerate component. Experimentally 
additional structures are seen which are not reproduced 
by the present calculations. It would therefore be inter- 
esting to include corrections to Migdal's theorem beyond 
the present calculations, and effects of the Coulomb in- 
teraction. 

This work has been supported by the Max-Planck- 
Forschungspreis . 
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TABLE I. Phonon level shifts and broadening from nu- 
merical estimates and Allen's formula. The numerical esti- 
mates are from fitting of q=0 spectral function (see text) 
and the estimates from Allen's formula use Eqs. (Ilq-hffl with 
A(0) = 6.2 per eV and g = 0.03 eV. 



V. APPENDIX: COEFFICIENTS c^f^ FOR 
RAMAN SCATTERING 

In the Raman scattering experiment (q = 0), the in- 
coming Hght with energy uji is scattered by losing the 
enerffii. of Uph to the phonon in sohd. Following Quang 
et al.t3, we treat the photon classically. The scattering 
rate Ip{uJi,LUph,) with a polarization vector p can then be 
expressed as 

Ip{uJi,UJph) OC y^^lmDa,a'iuJph)hPa{wi,UJph)h^,{Wi,UJph)*, 



(22) 



where h^{wi,ujph) is given by, as depicted in Fig. ||, 



n.rn.l 



X I dwG^(c^)G°„(c^ + iUpH)G'l{Lo + uui), (23) 

with matrix elements P^i of the dipole moment opera- 
tor at the polarization p between non-interacting elec- 
tronic eigenstates n and m. We consider the transitions 
iiu — > ilg^ since these transitions have large dipole matrix 
elements. Since a typical frequency of the light source is 
a few eV, and the ii„ — tig splitting is only about A^ ~ 1 
eV, transitions to higher states may, however, also play 
a role. 
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FIG. 3. Effective photon-phonon interaction vertex. The 
incoming photon (dashed lines) with frequency Wi scatters 
to a photon with frequency Wi — ujph by losing its energy to 
the a-th phonon mode (wavy line) with energy uiph- The 
interaction is mediated by the electrons (thick lines labelled 
as n, m, I) in the partially filled tiu band and in the empty 
tig band in the solid. The filled dot represents the elec- 
tron-phonon interaction with strength g",^ and the empty 
dots the electron-photon dipole interaction. 



with the bandwidth W ~ 0.5 eV, we ignore the Wi,ujph 
dependency in hP^{wi,ujph). Now factoring out the inte- 
gral, we have 



h^.{w.,u;p,)^l J2 Vjii^P^niPi 



i.ml ' 



(24) 



where the matrix elements P^i in Eq. (23) are reduced 

to matrix elements between local orbitals n, I at site i, 
pp 

We denote the three ii„ and tig-molecular orbitals as 
x, y, z, as in Ref. ^. For the first type of orientation, we 
can easily understand from symmetry of Ceo molecules 

that 






Pz 



pp 

i,yz 



P-- ^Lx=Py (orient. 1) (25) 



for a polarization direction p = {px,Py,Pz) with |p| — 1. 
Similarly for the second type of orientation, we only need 
to rotate the molecule 90 degree along z-axis, which gives 

Pr,.y = Pz, Kyz = -Py, Plz. = P. (orient. 2). (26) 

Note that -Pf„„ = 0. Therefore, leaving out common 
factors, we can write using Eq. (^, 

Ki ("^i ,^ph)=pl+pl- 2pl , pI+pI- 2pl 



tiP^{wi,ujph) = V3ipl-pl), -\/3(: 

}l%{Wi,UJph) ^iVipxPy, -2VipxPy 

hP^iwi.ujph) = '2VipyPzj 2VipzPx 

Kbi'^i^^ph) ^'^VipzPx, -2V3pyPz, 



Py-Px 



(27) 
(28) 
(29) 
(30) 
(31) 



where the first expressions on the right hand side are for 
molecules of orientation of type 1 and the second for type 
2. 

Finally, we take average of h^^UF- , over the polariza- 
tion directions, i.e., < h^^h^^, >= J ^h^^hP^,. Defin- 
ing c^, = I < h^ hF- , >, we obtain c as given in Eqs. 
(|I^]l|). Then 

< Ip{uji,u!ph) >pOc y^lmDa,aiujph), (32) 

a 

with Da.a defined in Eg. (||). 

Now, we look for the Dyson equation corresponding to 
the reduced Greens function D^^i,{q,uj) at q=0. we can 
rewrite Di^^ji, by splitting the coefficient matrix d^^, as 



We make an approximation to h^(wi,U!ph) where the 
intermolecular hoppings are ignored due to their rela- 
tively small bandwidth. Then, the indices in the inte- 
grand in Eq. (|2^ ) refer to orbitals on the same site. 
Since the three orbitals in ti„ or tig-orbitals are all 
equivalent and the integral depends on the indices only 
through their energies, the integral is only a function of 
Wi^ujpii, At. Furthermore, since Wi ^ W and ujph < W 



where r' is the 5x5 matrix defined as the unit matrix 



for o(i) = 1 and 



/ 1 \ 
-1 ' 
0-100 
0-1 

Vo 1 y 



for o{i) = — 1, 



^i\-l 



and r' 

the orbital indices, 



Now employing the matrix notation for 



k 



(34) 
(35) 
(36) 



with the reduced self-energy nipjv('jj) defined as 

'^Hii.jiy — / , Cm^,ila^' jV- (37) 

Finally by taking the summation over site indices i, j, we 
obtain the reduced phonon self-energy n^y(q = 0,w) for 
D^,y{q_^ 0,uj) as Eq. 
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